



MOTION OF A CIRCULAR CYLINDER AND n POINT 
VORTICES IN A PERFECT FLUID 



The paper studies the system of a rigid body interacting dynamically with point vortices in a perfect fluid. For arbitrary 
value of vortex strengths and circulation around the cylinder the system is shown to be Hamiltonian (the corresponding 
Poisson bracket structure is rather complicated). We also reduced the number of degrees of freedom of the system by 
two using the reduction by symmetry technique and performed a thorough qualitative analysis of the integrable system 
of a cylinder interacting with one vortex. 



In this paper we consider the system of a rigid cylinder interacting with point vortices. We start 
by indicating some known results on the subject from the classical hydrodynamics. For the most part, 
these results are presented in [18, 6, 9]. 

As far as we can see, Kirchhgoff [5] was the first who studied the dynamics of point vortices on a 
systematic basis. In particular, he obtained the equations of motion in Hamiltonian form and indicated 
integrals of motion. It was enough to show that the equations of motion governing the system of two 
or three vortices are integrable (the problem of three vortices was first analytically solved by Grobli) . 

The classics also considered the system of point vortices moving externally to rigid, stationary 
boundaries (the impermeable condition on the boundaries was assumed). Greenhill [12] studied in 
detail the motion of two vortices in a circular region. Havelock [13] investigated stability of n-gon 
stationary configurations in the region exterior to a circle. We should also mention Foppl's analysis 
of interaction of a rigid body with an ambient flow at low Reynolds numbers (more exactly, he 
investigated stability of the system of two vortices interacting with a circular cylinder embedded in a 
uniform flow). 

At the same time, the study of the dynamics of interacting vortices-body systems has been also 
a subject of interest in the classical hydrodynamics. It is known from the phenomenological theory 
developed by Prandtl [7] and Jukowski that when a body moves in a fluid, the thin boundary layer 
peels off the body thereby generating bound vortices. In their turn, these vortices exert a lifting force 
on the body which can be observed in aero- and hydrodynamical experiments. 

The problem of interaction of a rigid body and point vortices in a perfect fluid (the impermeable 
condition on the body's surface is assumed) can be studied within the framework of the Hamiltonian 
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mechanics. Various forms of the equations of motion for a rigid circular cylinder interacting with 
n point vortices have been recently (and practically simultaneously) obtained in [14, 15, 16]. The 
integrability of the equations in the case of one vortex (n = 1) was established in [10]. 

In this paper we will study in greater detail this case of integrability. We will also consider the 
simplest chaotic system of a cylinder and two vortices to which a reduction procedure will be applied 
resulting in reduction of degrees of freedom by one. 

Hereinafter, as in [10], the term rigid body will refer to a two-dimensional circular region. It 
should be noted that even in the case of an elliptic region the equations of motion become much more 
complicated and cannot be written in such a compact form. 



1. Hamiltonian form of the equations of motion 

The equations of motion for a cylinder and vortices with respect to a fixed coordinate frame Oxy can 
be written as [14] 

h = ~v+ grad ipi | r=r . , r c = v 

n ■ n ■ (1-1) 

avi = Xv 2 -'^2Xi(yi-yi), av 2 = -At>i+^ \j(xj-Xi), 

i=l i=l 

where r c is the radius- vector from O to the center of mass of the cylinder, v is the velocity of the 

cylinder, r\ is the vector from the center of the cylinder to the i-th vortex and = -^-n is the 

n 

vector from the center of the cylinder to the i-th inverse point (Fig. 1). Here R denotes the radius of 
the cylinder, the constant coefficient a involves the added mass of the cylinder; and the constants A 
and Aj are connected with the circulation around the cylinder and the vortex strengths by the formulae 

A = tr- , Aj = ^— . The density of the fluid is 2ir. 




Fig. 1 

The function <fi(r) represents that portion of the velocity potential <prr) which does not have a 
singularity at the point r = ri. The velocity potential in the region exterior to the cylinder reads 

2 n ~ 

^ = -^( r ^)- Aarct sI + Z] Ai ( arct g(f3| i ) " arct g(irr|))- (1-2) 

i=i 1 

Equations (1.1) were derived using a balance of linear momentum for the fluid within a circular 
boundary that encloses the body and the vortices. The fluid is assumed to be at rest at infinity [14, 16]. 

Thus the analysis of the system of a cylinder and vortices in a perfect fluid can be reduced to 
the analysis of a finite set of ordinary differential equations. It is easy to check that equations (1.1) 
preserve invariant measure, i.e., the divergence of (1.1) is zero. 
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As for the system of n point vortices [5], equations (1.1) can be shown to be Hamiltonian. 
Proposition 1. Equations (1.1) can be represented in the form 

Q = {Q,H} = Y,{^Ck}^- (1-3) 

k 



where Q are the components of the phase vector £ = (x c , y c , v 1 , v 2 , x 1 , y 1 , . . . , x n , y n ) and H is the 
Hamiltonian. For the components J^iC) = {Cn Cj} °f the structural tensor of the Poisson bracket 
structure the Jacobi identity holds: 

( dJ jk dJ in - dJk\ 
\ [ Ja ^Q + Jkl ^ + Jjl ~d^ ) = ° Vi ' j ' k 

Proof. 

Equations (1.1) has an integral of motion: 



H = \ av 2 + \ J2 (a? ln(rf - R 2 ) - \\ In rf) + \ ^ \\ 3 In ^^'^ + - - j 



(1.4) 



'J 



This integral resembles the Hamiltonian of the system of n point vortices [5]. 

Let H be our Hamiltonian. Now we have to find a skew-symmetric tensor J- sue 
equations of motion (1.3) coincide with (1.1). The non-zero components of this tensor are 



s , 1 4-R 2 (x 2 -y 2 ) l2R 2 x m 
{vi, Xij = — , {vi, yij 



r 4 ' «-"x, a 4 

i i 



1 2R 2 x m 1 rf + R 2 (x 2 -y 2 ) 

{V 2 , Xij = — , {^2, yij ~ 



a rf — a rf (1.5) 

r t A ^i T i ~ [ 1 1 

{vi,v 2 } = ^-)^ 2 — , {x i ,y i } = — r , 



a" a rj *i 
{x c , vi} = {y c , v 2 } = a" 1 . 



It is easy to verify that for the Poisson bracket (1.5) the Jacobi identity holds. ■ 

The Poisson bracket structure (1.5) is non-degenerate, therefore, by the Darboux theorem, it can 
be reduced to a canonical form ({<^, pA = <5y). However, in our further analysis canonical coordinates 
will not be used. 

The Lie-Poisson bracket structure for (1.1), (1-2) under the condition A = — was studied 
in [16]. In this work stability of an equilibrium configuration of the system of two vortices behind a 
steadily moving cylinder is discussed. The Poissom bracket (1.5) was first obtained in [10]. 

The fact that the general equations are Hamiltonian is not a priori obvious and does not seem to 
follow from the Lagrangian formalism. The Hamiltonian form of the equations allows us to apply the 
highly developed perturbation theory techniques (e.g. the KAM theory) and other specific methods 
of qualitative analysis. The Liouville theorem on integr ability and its geometrical extension suggested 
by Arnold [1] can be applied to the analysis of these equations. 

2. Problem of advection 

The problem of finding pathlines of the fluid for a given motion of the cylinder r c = r c (t) and the 
vortices r { = Vi(t), i = 1, . . . , n is known as the problem of advection. As mentioned above, the 
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ambient flow is potential with potential (1.2), therefore the equation of motion for a passive particle 
with respect to the cylinder-fixed frame of reference looks like 

r = grad<p(r)| (t) , r = (x, y). (2.1) 

Obviously, equations (2.1) are Hamiltonianian with respect to the standard Poisson bracket struc- 
ture {x; y} = 1; the Hamiltonian is time-dependent and coincides with the stream function for the 
potential (1.2), that is, 

N (2.2) 
+ M r - ri(t)\ 2 -\n\r-n(t)\ 2 ) . 

i=l 

3. Symmetry and integrals of motion 

The equations of motion (1.1) are invariant under the action of the Euclidean group E(2), therefore, 
by Noether's theorem for Hamiltonian systems, there exist three integrals of motion. 
The integrals 

Q = av 2 + \x c - \{Xi ~ Xi), P = aVl - Xy c + ^ X^ - Vi ) (3.1) 

correspond to translations along the coordinate axes. These integrals are a generalization to the 
classical linear momentum. On the other hand, the vector (Q, P) can be considered as a counterpart 
of the vector of the center of vorticity of the system of n vortices and cylinder. For the system of n 
vortices this notion was introduced in [5, 6]. For A / 0, the origin of coordinates can be so chosen 
as to make Q = P = 0, meaning that the center of vorticity can be always shifted to the origin of 
coordinates. 

The third integral, corresponding to the invariance under rotations about an axis perpendicular 
to the plane of motion, is 

n n / \ 

/ = a(v lVc - v 2 x c ) - \\rl - \ \rf + \ £ \\ - 1 (r <9 r c ). (3.2) 



2 ^ 1 1 2 

i=l i=l 



In the paper [10], an additional integral for equations (1.1) but with the dynamics for the center 
of the cylinder excluded was indicated. This integral looks like 



-y i v 1 )+{X i -X)r 2 +X i ^ ) 



F = a 2 v 2 + A * ( 2a ( 1_ §) ( Xi 

+2^A,A,(n,r-,)(l-f)(l-f) 



(3.3) 



i<j v » y v 3 

The integrals F,I, P and Q satisfy the equation 



n 

2 



F = 2X1 + P 2 + Q 2 + 2R 2 Y A 



i=i 

Comment. The Hamiltonian vector field that corresponds to the integral (3.2) looks like 

x i = {C 1} = (j/o ~x c , v 2 , -v 1: y lt -x lt . . . , y n , -x n ) (3.4) 
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The Poisson bracket of the integrals Q, P, and / differs from the Lie-Poisson bracket for the e(2) 
algebra by a constant (co-cycle [1]), that is, 

{Q,P} = \, {I,Q} = P, {I,P} = -Q. (3.5) 

Thus, for A 7^ the number of degrees of freedom of the system governed by (1.1) can be reduced 
by two and even by three if A = and P = Q = 0. 

Corollary 1. The dynamics of a cylinder and one vortex is integrable in the Liouville sense. 

Corollary 2. For A = and P = Q = 0, the dynamics of a cylinder and two vortices is integrable 
in the Liouville sense. 

4. Complex form of equations of motion and the Dirac bracket 

As mentioned above, in [10] reduced equations were used, i.e., the equations we dealt with were exactly 
equations (1.1) but without the equation r c = v. The elimination of this equation from the general 
system (1.1) now can be interpreted as the reduction by symmetry due to the integrals Q and P (3.1). 
For A 7^ the reduction can be carried out by restricting the dynamics to a joint level surface of the 
integrals (3.1). Since A 7^ 0, we assume that the origin of the fixed coordinate frame is at the center 
of vorticity, i.e., P = Q = 0. Then, we substitute r c for v i in (3.1). The first order equations in the 
positional variables result. It is convenient to write these equations in the complex form 



av 



i\z c + i x j 



3=1 



| fe = _ c + ^ +i A _ a ^ + V A _1 1\ (4.1) 

z k = -=— k = 1, . . . , re, 

z k 

where z c = x c + iy c and z k = x k + iy k define the position of the cylinder's center and the vortices, and 
v = v 1 + iv 2 is the velocity of the cylinder's center. 

Obviously, equations (4.1) are Hamiltonian. The Hamiltonian can be obtained by replacing the 
cylinder's velocity in (1.4) with the expression in the right-hand side of the first equation of (4.1) 

H = 2a" ( Xx ~ 12 Xi ^i ~ x ^> ) + ^[ x y -^2 x i(yt-yi)j + 

+| E Of Hri ~ R 2 ) ~ A,A In rf) + \ ^ \\j In - - ^ ^ + ^ . 

i i<j ' 1 3 * 

The Poisson bracket for (4.1) is 

{x c ,y c }=y {Xi,yi} = -j-, i = l,...,n. (4.3) 

This bracket can be obtained via the Dirac reduction procedure [11], which consists in restricting the 
bracket (1.5) to the manifold Q = P = 0. In other words, for A 7^ 0, the cylinder can be considered as 
an (n + l)-th "compound vortex", and equations (4.1) govern the system of re + 1 vortices. 
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Comment. The Dirac bracket {-,-}d on a manifold Af c which is a level surface fi(C) = c i-> °i = const, 
i = 1, . . . , k is given by the formula 

{g, h} D = {g, h} + ^{.9, fi}cij{h, fj}, 

where ||cy|| = /j}!! -1 and {•, •} is the original Poisson bracket. In our case we have 

{g, h} D = {g, h}-± ({g, Q}{h, P} - {g, P}{h, Q}) , 

where {•, •} is the Poisson bracket (1.5) 

Equations (4.1) are invariant under rotations about the origin (the center of vorticity) and there- 
fore have an additional integral of motion: 



t=i 



(4.4) 




-62,7 



-30.1 



Fig. 2. Motion of a cylinder (solid line) 
and a vortex couple (ri = — T2); a = 

= 20, r = 8, r x = 10, r 2 = -10, r = i, 

F = 2; the initial conditions (relative to 
the center of the cylinder) are i>i(0) = 0, 
v 2 (0) = 0, xx(0) = -3, Vl (0) = 0, x 2 (0) = 
= 3, ife(0)=0. 



This integral can be obtained upon substitution of the expres- 
sion for r c from (4.1) for v in the integral (3.3). 

Now we can formulate some properties of the motion of 
the system of a cylinder and n vortices. These properties are 
analogous to those formulated by Synge [17] for the system of 
n vortices. 

Proposition 2. Suppose that the vortex strengths are of 
the same sign and AAj < 0; then the trajectory of the cylinder's 
center and the vortices belong to a bounded region. 

Proposition 3. Suppose that A 7^ and the vortex mo- 
tions relative to center of the cylinder are bounded (i.e., \/i 
I (t) I is a bounded function), then the absolute motion of the 
cylinder is also bounded. 

Here and in the sequel, a motion with respect to the 
fixed frame will be referred to as an absolute motion and with 
respect to the center of the cylinder as a relative motion. 



Comment. It should be noted that for n = 1 and A 7^ a more 
general statement is valid: the absolute motion of the cylinder is 
bounded if and only if the relative motion of the vortex is bounded. 
Simulations have shown that already for n = 2 this is not true: the 
trajectories of the cylinder and two vortices for the case Ai = — A2, 
|A| < I A x I are shown in Fig. 2. As might be expected, eventually the vortices start drifting to infinity while the 
cylinder moves along a curve which gradually takes the shape of a circle. 



5. Motion of a cylinder and one vortex 

Integrability of the equations of motion. Reduction to a system with one degree of 
freedom. Consider in greater detail the system of a cylinder interacting dynamically with only one 
vortex, i.e., n = 1. Let r\ = r = (x,y). With the assistance of the first integrals (3.1) and (3.2), the 
solution to the equations of motion follows in terms of quadratures. Using the integrals, our system 
can be reduced to a system with one degree of freedom. The algebraic reduction that we will now use 
is analogous to the Routh reduction. 
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As the variables of the reduced system, it is reasonable to choose integrals of the field of symmetry 
vi (3.4) (see, for example, [2]). We put 

pi = a(xvi + yv 2 ), p 2 = a(xv 2 - yvi), p = x 2 + y 2 . (5.1) 
The Poisson brackets for these variables are 



{pi,P2} = (A - Ai)p + f (pi + (p 2 - XiR 2 ) 2 ), 

Ai 

{ Pl ,p} = 2p + ^{p 2 - XiR 2 ), {p 2 ,p} = -2^-. 



(5.2) 



In terms of (5.1) the integrals (1.4) and (3.3) read 



H = + \\\ In (p - R 2 ) - iAxAlnp, (5.3) 

F = *1±& + 2Ax ( 1 - — ) p 2 + A? (p + — ) - AA 1P . (5.4) 



P V P J VP 

The rank of the Poisson bracket structure (5.2) is two, i.e., the structure is degenerate. Hence, 
the reduced system has one degree of freedom. One variable of the set (p±, p 2 , p) can be eliminated 
using the integral (5.4), which is a Casimir function for the structure (5.2). 

Traditionally [1], on a two-dimensional level surface of the Casimir function (a symplectic leaf), 
canonical coordinates ({q, p} = 1) are introduced. The reduced equations are 

■ dH . dH 

dp ' dq 

In contrast to this traditional approach, the local coordinates that we will use account well for 
the leaf's geometry but are not canonical. The phase portraits in terms of these coordinates are vivid 
and illustrative. Our attempts to find a simple and natural set of canonical coordinates have not been 
successful. 

General properties of motion. The reduced system (5.2) describes the motion of the vortices 
relative to the center of the cylinder. Before proceeding to a qualitative analysis of the reduced system, 
we will indicate some general properties of the absolute motion. 

Proposition 4. The absolute motion of the cylinder is bounded except maybe for the two cases: 

1. A = Ai; 

2. A = 0; 

Comment. Condition (1) corresponds to the case where the circulation around the cylinder is zero, and 
condition (2) represents the case where the circulation around the cylinder and the vortex is zero. 

Proof. 

Assume the contrary: A ^ Ai, A ^ and the motion of the cylinder is not bounded. Then, 
according to Prepositions 2 and 3, p is an unbounded function of time and AAi > 0. Dividing (5.4) 
by p yields 



Obviously, 

, , 2 / Tt-i \ 2 

lim 

p— >oo 



AAi. 
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Therefore, there exist functions a(p) = o(l) and (3{p) such that Ai + = \J\\X + a cos /?, = 

= \/AiA + asin/3. 

Upon substitution of these expressions into (5.3), we get 

2i J= f fc + A 1 ln("-f )A ' 
za 

The second term in the right-hand side is of 0(lnp), the factor k is separated from zero, that is, 
k = (Ai — V AiA) 2 + a + 2Ai(\/AiA — \A\iA + a cos /?) > k\ = const > 0. This is a contradiction with 
the fact that H is constatnt. ■ 

Comment. One can easily note that for A = Ai and A = an unbounded motion of the cylinder always 
exists. Therefore the above preposition serves as a criterion of boundedness of the absolute motion. Most likely 
this criterion remains valid for the case of an arbitrary (finite) number of vortices, but this is not proved yet. 

The statement given below provides a nice " geometro-dynamicai" insight into the structure of 
absolute motions of the cylinder. 

Theorem. Suppose that A / 0; then for each periodic solution of the reduced system (5.1) there 
exists a rotating coordinate frame with the origin at the center of vorticity such that with respect to 
this frame the vortices and the cylinder move along closed curves. 

Proof. 

It follows from (4.1) that 

yxi - xy x = —, yyi + xx x = -— + —(R - p). (5.6) 

By the assumption of the theorem, the functions in the right-hand side of these equations are periodic 
functions with period T. Therefore 

;:) ♦•(;)• *<X%) 

where $i and $2 are periodic functions with period T. Substituting these expressions into (4.1), we 
get 

ax = \y — yG2 + xGi, ay = —\x + yd + XG2, (5-8) 

where G± and G2 are also T-periodic functions. Let A be the matrix for the system of linear differen- 

t t 

tial equations (5.8). It can be easily verified that A • J Adt = J A dt ■ A, hence (see, for example, [4]) 



the fundamental matrix, X, of equations (5.8) can be written as 

t 

X = eo 



Let us represent A as a sum 



A = B + C, 



B = l( A - <G 2 ) \ c= i( Gi -G 2 + (G 2 ) 



a \ -A + (G 2 ) J ' ~ a \G 2 - (G 2 ) G 1 

Note that (Gi) = 0, otherwise equations (5.8) have unbounded solutions which is in a contradiction 
with Preposition 3. Since the matrices B and C commute, we have 

A — (G2) , . A— (Gg) 
cos i sin i 

■ A-(G 2 ), A-(C7 2 ) + 
sm 1 cos 1 



456 



REGULAR AND CHAOTIC DYNAMICS, V. 8, JVM, 2003 



MOTION OF A CIRCULAR CYLINDER AND n POINT VORTICES IN A PERFECT FLUID (g£S 



Here G is a T-periodic matrix. Thus, with respect to a coordinate frame (with origin at the center of 
vorticity) rotating at a rate (A — {G2))/a the cylinder moves along a closed curve. Let us prove now 
that the vortex also moves along a closed curve. Let x*,y* be the coordinates of the vortex in the 
fixed coordinate frame. Then x* = x c + x, y* = y c + y. It follows from (5.7) that 

x :)=(* + e).( x 

y J V v 

A - (G 2 ) . A - (G 2 ) \ 
cos i^-t sin \ , / x 0) 

" sm a t cos a */ 

Here E is the identity matrix, and the matrix (<fr + E)G is T-periodic. 




Fig. 3. a) Trajectories of the cylinder (solid line) and the vortex (here the two frequencies of the motion are 
rationally related), a = 6, T = —3, I\ = 1, R = 1, F = 2, ; the initial conditions (relative to the center of 
the cylinder) are v%(0) = 0.24, Va(0) = 0, x(0) = 0, y(0) = 1.31; b) The trajectories shown in diagram a in the 
frame of reference rotating at a rate 0.029 (the origin, + , is at the center of vorticity). 

An absolute motion of the cylinder (solid line) and the vortex (dashed line) is shown in Fig. 3 a. 
In Fig. 3 b this motion is shown in a rotating frame of reference whose origin, 0+, is at the center of 
vorticity. The value of the physical parameters for this motion are a = 6, T = — 3,Ti = 1, R = 1, 
F = 2 and the initial conditions are «i(0) = 0.24, u 2 (0) = 0, x(0) = 0, y(0) = 1.31. 

Qualitative analysis of the reduced system. The geometry of the symplectic leaf (5.4) of 
the Poisson bracket structure (5.2) (the phase space of the reduced system) is governed by AAi. The 
symplectic leaf is compact if AAi < and non-compact if either AAi > or A = 0. Let us consider 
these three cases. 

1. Compact case (AAi < 0). From (5.4) it follows that 

{\ 1 p + p 2 -\ 1 R 2 ) + p\ + (^\ lp + K) =K 2 : K= ^7===- (5.9) 

^ ' 2y— AAi 

The symplectic leaf is diffeomorphic to a two-dimensional sphere. For real motions F ^ 2Aiii 2 (Ai — A). 
Since the phase space is compact, p is a bounded function, meaning that the distance between the 
vortex and the cylinder cannot grow infinitely. Moreover, in view of Preposition 3, the absolute motion 
of the cylinder also is bounded. 



REGULAR AND CHAOTIC DYNAMICS, V. 8, J\M, 2003 



457 



A.V.BORISOV, I. S. MAMAEV, S. M. RAMODANOV 



On the leaf we intoduce coordinates p and ip by the formulae 

sin ip cos ip, 

(5.10) 



X\p + P2 — X±R 2 = K cos (p cos ip, p\ = K sin <p cos ip, 
py/^XX 1 + K = Ksmip, ipe (-§,§) , <P& [-7r,7r] 



Comment. For a fixed value of K, the coordinate %h satisfies the relation p— ■ — > R 2 . 

To find stationary solutions of the reduced system, consider the differential equations in pi,P2 
and p: 



Pi = {pi,H} 



XiR^Xap - 2R 2 X 1 Xap 2 - X\ap 3 + XlR 2 ap 2 + AiAap 3 



+ 



p 2 a(-p + R 2 ) 

(-R 4 Xip + Xap 2 + - XiR 2 p 2 - XR 2 ap + XR 2 p 2 - X iap 2 + Aip 3 - Xp 3 )p 2 

p 2 a(-p + R 2 ) 

(p 2 -R i )p 2 2 + (p 2 -2R 2 p + R*)p 2 

+ f? a (-p + R2) 

(Xp 3 -X 1 p 3 +XR 2 ap-R 6 X 1 +X 1 ap 2 -Xap 2 -XR 2 p 2 +X 1 R 2 p 2 +R' l X 1 p)p 1 

p 2 a(-p + R 2 ) 

{-2R 2 p + 2R A )p 2 p 1 



(5.11) 



P2={P2, H}= + 



+ 



p 2 a( _ p + jR2) 

(-2p + 2i? 2 )pi 



a/3 



The last equation implies pi = and the right-hand side of the second equation becomes zero. There- 
fore, on the phase portrait, all stationary solutions belong to the lines ip = p = and p = ±7r. 

To determine p 2 and p we should use (5.9) and equate to zero the right-hand side of the first equation. 
The number of stationary solutions depends on F. 

Equations (5.11) remain unaltered under the transformation (p\, p 2 , p, t, A,Ai) — > (—pi, — p 2 , 
p, —t, — A, — Ai). In view of (5.10), this transformation implies the following change of variables and 
physical parameters (p, ip, X, X\) — ► {p + ir, ip,—X,—X\). Therefore, a change in sign of the circulation 
and the vortex strength results in a shift of the phase curves along the <p-axis by ir. Thus, with no 
loss of generality, we can assume that A < and Ai > 0. 

Phase portraits for various values of F ( the value of the other parameters are fixed a = 20, 
r = —1, Ti = 0.5 and R = 1) are shown in Figs. 4 a, 4 b and 4 c. The "non-physical" area (p R 2 ) is 
shown in grey. 

The solution curves of the reduced equations are the level curves of (5.3) on the surface (5.4). 
The level curves are closed, hence the solutions of the reduced equations are periodic. 

The fixed points on the phase portraits (stationary solutions) represent a motion in which the 
cylinder and the vortex move along concentric circles whose centers are at the center of vorticity. For 
example, the trajectories of the vortex and the cylinder that correspond to the elliptic point in Fig. 4 
are shown in Fig. 5. With an increase in F two more fixed points appear, the motion corresponding 
to these points is qualitatively shown in Fig. 6. 

2. Non-compact case (AAi > 0). It follows from (5.5) that 

2 2X 2 R 2 F 

{X lP + p 2 -X 1 R 2 ) 2 +p 2 - (^/XXlp-K} =-K 2 , K= 1— . (5.12) 
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0,5tt * i> 




Fig. 4. Phase portraits of the reduced system in the compact case (AAi < 0): a = 20, r = — 1, Ti = 0, 5, R = 1. 





0> 



Fig. 5 



Fig. 6 



Symplectic leaves are diffeomorphic to a hyperboloid of two sheets. For real motions F>2\iR 2 (Xi— A). 
On symplectic leaves we introduce local coordinates: 

Xip + p2 — X\R 2 = K cosh (^sinh-^, pi = K sinhip, 
p\J AAi — K = —K cosh tp cosh y?, 



(5.13) 



Arguing as in the compact case, we can assume that A > and Ai > 0. In view of (5.11) and (5.13), 
on the phase portrait, all fixed points lie on the axis (p = 0. It is interesting to note that, unlike the 
compact case, the topology of the phase portrait is determined by the sign of the difference A — Ai and 
does not change qualitatively as the constant F varies. Typical phase portraits are given in Fig. 7 a 
(A = Ai), 7 b (A < Ai) and 7 c (A > Ai). 

The phase curves are closed, hence the solutions of the reduced equations are periodic functions 
of time. The only exception is the case A = Ai (Fig. 7). As shown above, only in this case the cylinder 
may have an unbounded motion (Preposition 4). 



REGULAR AND CHAOTIC DYNAMICS, V. 8, J\M, 2003 



459 




-9,3 1 

r = 12 



Fig. 7. Phase portraits of the reduced system in the non-compact case (AAi > 0): F = 25, a = 20, Ti = 10, 
R= 1. 

3. Non-compact case (A = 0). In the system of equations (1.1), the equations governing the 
motion of the vortex are uncoupled from the first equation. Indeed, it follows from (3.1) that 

avi = -Aiyi/(xi, yi) + ci, av 2 = XiXif(xt, yi) + c 2 , 

where f{x\,y\) = — l + R 2 /{x\ + y\) and c\, c 2 are arbitrary constants. Substituting these expressions 
into the first equation (1.1), we obtain a closed system of equations in the unknowns xi and y\. The 
solution curves of this system coincide with level curves of the integral (5.3). It can be shown that the 
level curves are closed, hence x\ and y\ are periodic functions of time. The evolution of the center of 
the cylinder is governed by the equations 

x = a (~ Ai (yi/( x i'yi)) + c i)t + 9i(t), 

y = -(\ 1 (x 1 f(x 1 ,y 1 )) + c 2 )t + g 2 {t), 

where gi(t), g 2 (t) are periodic functions. Thus, there exists a uniformly moving coordinate system in 
which the orbits of the cylinder and the vortex are closed. 
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6. The case of two vortices 

Suppose now that n = 2, i.e., we are going to consider the system of a cylinder interacting dynamically 
with two vortices. By analogy with the case n = 1, we use the integral (3.3) to reduce the number 
of degrees of freedom by one. As before, we take some quantities invariant with respect to rotations 
about the cylinder's center as the new variables, namely, 

pi = a(xivi + yiv 2 ), P2 = a(xiv 2 - yivi), Pz = x x x 2 + yiy 2 , P4 = xiy 2 - x 2 yi, 

2,2 2,2 

n=x 1 + y 1 , r 2 = x 2 + y 2 . 
The Poisson brackets for these variables are as follows: 



{Pl,P2> 
{Pl,Ps} 
{Pl,Pi} 



(pi - 2X 1 R 2 p 2 r 2 + (i? 4 - r\)\\r\ + r\\\\r\ + r^X^R 4 - r 2 ) 



\\r\r 2 

rl(P3P2 ~ PaPi) + Ai(-(i? 2 - ri)rfp 3 + 2R 2 p\r x - R 2 r\r 2 + r\r\) 

\\r\r 2 

r 2 (piP3 +P2P4) - \iPi{2R 2 p2,n + {R 2 - ri)r|) 



1 2 

?2 



2p 2 - 2(R 2 - n)Ai r . 2p 3 (R 2 -r 2 ) 



{pi,n}= r » {Pi' r 2} 

{P2,P3> 



r 2 (PiP3 + P2P4) ~ \m(2R 2 p 3 ri - (R 2 + r{)r%) 



{P2,Pi} 



\ 2 ' 

\\r\r 2 

r 2 (p3P2 ~ P4Pi) + Ai(-(i? 2 + ri)r|p 3 - 2i2 2 p|ri + ,R 2 rfr2 + rfr|) 



'2 

r , 2pi 2p 4 (i? 2 - r 2 ) \ 2 r 2 - \ lTl 

{P2,n} = - — , {P2,^2} = ~ 2 , {P3,P4j = , 

{P3,ri} = -^, {^3,^2} = -^, 

{P4,n} = — 7—, {P4,r 2 } = ^— , |ri,r 2 } = 0. 

Ai A 2 

This Poisson bracket structure is degenerate, and its rank is four. The integral (3.3) and the 
obvious relation p 2 +p 2 = T\r 2 are Casimir functions for this structure. For the reduced system to be 
integrable, one more first integral is needed. 

To explore the reduced system numerically, we have used the Poincare surface-of-section technique. 
The variables r 2 , P3, P4 can be considered local coordinates on a three-dimensional manifold on which 
the integrals (1.4) and (3.3) are constant and the equation p\+p\ = r±r 2 is fulfilled; p^ is the cross 
variable. Two Poincare surface-of-section plots are given in Fig. ??. The chaotic behavior of solutions 
proves that in the general case an additional first integral does not exist. 



7. Conclusion 

Very often equations of motion that have not been derived within the framework of Lagrangian for- 
malism (i.e., using the calculus of variations) are not Hamiltonian (even though the energy may be a 
conserved quantity for these equations). For example, for equations of the nonholonomic mechanics 
there are dynamical obstacles preventing the existence of a Poisson bracket structure [3]. Since the 
equations of motion (1.1) are Hamiltonian, they have the generic features of Hamiltonian dynamical 
systems: for any value of parameters there are no attractors (e.g., strange attractors) in the phase 
space; at the same time, there are invariant KAM tori separated with stochastic layers. In this paper, 
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a) 



b) 



Fig. 8. Poincare surface-of-section plots for the system of a cylinder and two vortices: a) H = 7, F = 5, a = 9, 
r = 1, T 1 = T 2 = 8, R = 1, Pi = -5; 6) H = 10, F = 60, a = 4, T = 2, r x = -1, T 2 = 10, i? = 1, p 4 = 0. 

the chaotic system of a cylinder and two vortices has been touched briefly. It seems that it would 
be interesting to explore some particular motions of this system (both regular and chaotic) in greater 
detail. 

Mention should be made of the paper [8] where a modification of the famous Bjerknes problem, 
the system of two dynamically interacting 2D cylinders, is considered. The equations of motion for 
this system are not integrable. 

This work was supported in part by Leading Scientific School of Russia Support grant 
no. HIH-136.2003.1. 



References 



[1] V. I. Arnold. Mathematical methods in classical me- 
chanics. Moscow: Nauka. 1991. (In Russian) 

[2] A. V. Borisov, I. S. Mamaev. Poisson structures and 
Lie algebras in the Hamiltonian mechanics. Izhevsk, 
RCD. 1999. (In Russian) 

[3] A. V. Borisov, I. S. Mamaev. (Strange attractors in 
rattleback dynamics) UFN. V. 173. JVM. P. 407-418. 

[4] B. P. Demidovich. Lectures on the mathematical the- 
ory of stability. Moscow: Nauka. 1967. (In Russian) 

[5] G. Kirchhoff. Vorlesungen iiber mathematische 
Physik. Mechanik, Leipzig. 1874. 

[6] H. Lamb. Hydrodynamics, Ed. 6-th. N. Y. Dover publ. 
1945. 

[7] L. Prandtl. Foundamentals of Hydro and Aeromechan- 
ics. Dover Publ., Inc. New York. 1957. 

[8] S. M. Ramodanov. Motion of two circular cylinders 
in a perfect fluid. In: Fundamental and applied 
problems in the theory of vortices, eds. A. V. Borisov, 
I. S. Mamaev and M. A. Sokolovskij. Moscow-Izhevsk: 
Institute of Computer Science. 2003. P. 327-335. 

[9] P. G. Saffman. Vortex Dynamics. Camb. Univ. Press. 
1992. 

[10] A. V. Borisov, I. S. Mamaev. An integrability of the 
problem on motion of cylinder and vortex in the ideal 
fluid. Reg. & Chaot. Dyn. 2003. V. 8. P. 163-166. 



[11] P.A.M.Dirac. Generalizated Hamiltonian Dynam- 
ics. Canadian Journal of Math. 1950. V. 2. 
JYs2. P. 129-148. (See also [Borisov A.V., Ma- 
maev I. S. Dirac bracket in geometry and mechan- 
ics, in Dirac P. A. M. Lectures on theoretical physics. 
Moscow-Izhevsk: RCD, 2001, in Russian) 

[12] A. G. Greenhill. Plane vortex motion. Quart. J. Pure 
Appl. Math. 1877/78. V. 15. Ns58. P. 10-27. 

[13] T. H. Havelock. The stability of motion of rectilinear 
vortices in ring formation. Phil. Mc. 1931, Ser. 7. 
V. 11. JV»70. P. 617-633. 

[14] S. M. Ramodanov. Motion of a Circular Cylinder and 
a Vortex in an Ideal Fluid. Reg. & Chaot. Dyn. 2001. 
V. 6. J\VL. P. 33-38. 

[15] S. M. Ramodanov. Motion of a Circular Cylinder and 
N Point Vortices in a Perfect Fluid. Reg. & Chaot. 
Dyn. 2002. V. 7. m. P. 291-298. 

[16] B. N. Shashikanth, J. E. Marsden, J. W. Burdick, 
S. D. Kelly. The Hamiltonian structure of a 2D rigid 
circular cylinder interacting dynamically with N point 
vortices. Phys. of Fluids. 2002. V. 14. P. 1214-1227. 

[17] J. Synge. On the motion of three vortices. Can. J. 
Math. 1949. V. 1. P. 257-270. 

[18] H. Villat. Lecons sur la theorie des tourbillions. 
Gauthier-Villars. 1930. 



462 



REGULAR AND CHAOTIC DYNAMICS, V. 8, JVM, 2003 



